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TRANSIENT RESPONSE OF SHELLS OF REVOLUTION BY DIRECT 
INTEGRATION AND MODAL SUPERPOSITION METHODS 

By Wendell B. Stephens and Howard M. Adeiman 
Langley Research Center 


SUMMARY 

This report describes the results of an analytical effort to obtain and evaluate 
transient response data for a cylindrical and a conical shell by use of two different 
approaches: direct integration and modal superposition. The purpose of the investiga- 
tion was twofold: to provide reliable data for use in the checkout of general-purpose 
codes and to gain some insight into the relative advantages of each technique. 

Results are obtained for the cylindrical shell which is freely supported at both 
edges and for the conical shell which is freely supported at the small diameter and free 
at the large diameter. Both shells are subjected to an initial velocity excitation. Both 
the direct integration analysis and the modal superposition analysis included such second- 
ary effects as rotary inertia and transverse shear deformation. As a result, the impor- 
tance of these effects as well as the ease of their inclusion in the calculations were eval- 
uated. In addition, the numerical integration analysis included nonlinear effects due to 
large displacements. Consequently, the importance of these effects was evaluated also. 

The principal conclusions drawn from the present study are that for the types of 
problems investigated herein, the inclusion of nonlinear terms is more important than the 
inclusion of secondary effects although there are thin -shell structures where these 
secondary effects will be important. The advantages of the direct integration approach 
are that geometric nonlinear and secondary effects are easy to include and high-frequency 
response may be calculated. In comparison to the modal superposition technique the 
computer storage requirements are smaller. The advantages of the modal superposition 
approach are that the solution is independent of the previous time history and that once 
the modal data are obtained, the response for repeated cases may be efficiently computed. 
Also, any admissible set of initial conditions can be applied. 

INTRODUCTION 

Dynamic response behavior can be an important design consideration for thin shells 
of revolution which are often used as structural members in aerospace vehicles and civil 



structures. As a result, there have been attempts in recent years to provide analytical 
procedures for accurately calculating dynamic response of general shells of revolution. 
The two methods that have shown the most promise for a wide variety of shell geometries 
are direct integration of the equations of motion by using finite differences (refs. 1 to 5) 
and modal superposition based on a finite -element calculation of the modal character- 
istics of the shell (refs. 6 to 8). Each approach is most advantageous for a certain sub- 
class of the general problem. 

Recently, the trend in analytical methods development has been toward large 
general-purpose computer programs (refs. 9 and 10). One important task in the develop- 
ment of such programs is comparing those results with the results of special-purpose 
codes in order to generate confidence in the minds of potential users. Thus, the purpose 
of the present paper is to provide some reliable results for transient response of shells 
of revolution for use in checking more general-purpose computer programs. 

As a part of the present work, it was desirable to carry out some extensions to 
previous work and to develop some new procedures. These extensions are presented in 
three appendixes. Appendix A contains modifications to the direct integration procedure 
of reference 1 to include the effects of rotary inertia and transverse shear deformation. 
Appendix B describes the modal superposition algorithm. Appendix C describes the 
finite -element derivation for modal characteristics of cylindrical shells with rotary 
inertia and transverse shear effects included. A fourth appendix, appendix D, lists in 
tabular form all numerical results shown in the figures. This is provided so that mean- 
ingful checks with other analyses can be made. 


SYMBOLS 

j~A^j matrix which relates coefficients of assumed displacements to displacements 

and rotations at ends of element 


a Q k ,aj ^, a 2 k > a 3 k coefficients in assumed displacement function for u 

b 0 k’ b l k’ b 2 k’ b 3 k coefficients in assumed displacement function for w 

Eh 

C membrane stiffness, - 

1 - v z 

[C k ] matrix defined in equation (C21) 


C 0 k ,Cj ^jCg k ,Cg k coefficients in assumed displacement functions for /3 


2 



D 

E 


e l’ e 2 

W 

G 


h 

i = fi 
K 



k 


2 


N 

K] 


n 





r 


flexural stiffness, 
elastic modulus 
meridional and circumferential strains, respectively 
matrix defined in equation (C12) 

shear rigidity, ^TTv ) 

shell thickness; H in computer tables 

total number of elements 

stiffness matrix for nth harmonic 

shear correction factor (0.87 for a thin shell) 

meridional length 

element mass matrix 

mass matrix for the nth harmonic 

meridional moment resultant; M-X in computer tables; MX in computer plots 

membrane force resultants 

positive integer; N in computer tables 

transverse shear stress resultant 

normal coordinate of the system 

matrix defined in equation (C7) 

radial distance; R in computer plots and tables 


Eh 

12(l - i/ 2 ) 
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M 

s 

s k 

T 

t 

U,V,W 

u,v,w 

V 

v o 

[x] 

X 

M 

W 

N 



0 

M 


^13 


element stiffness matrix 

meridional coordinate; S in computer plots and tables 
distance from first edge of cylinder to center of kth element 
kinetic energy 
time; T in computer plots 

meridional, circumferential, and normal displacement components at any 
point in shell/ respectively 

meridional, circumferential, and normal displacement components of point 
on middle surface, respectively 

strain energy (appendix C) 

initial velocity 

matrix defined in equation (Cl 7) 
meridional coordinate in a finite element 
matrix defined in equation (C25) 

vector containing all unknown displacements and their spatial derivatives 
matrix defined in equation (C27) 

coordinate measured normal to middle surface of cylinder 

rotation of a point on middle surface of cylinder 
vector defined in equation (C16) 
transverse shear strain 
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length of kth element 


€ k 

W 


'n 


diagonal matrix of eigenvalues for nth harmonic (appendix B) 


eigenvalue of nth harmonic (appendix B) 


v Poisson’s ratio 

p mass density 

<y l > (J 2> (y X3 stress in the meridional, circumferential, and transverse directions, 

respectively 

matrix of eigenvectors for nth harmonic (appendix B) 

eigenvector for nth harmonic (appendix B) 

cone angle (fig. 1(b)) 

oj angular frequency, rad/sec 

A prime indicates a derivative with respect to the meridional coordinate s. 
A dot indicates a derivative with respect to time t. 


ANALYTICAL METHODS 


Direct Integration 

The results in this paper illustrating the direct integration method are obtained 
from the shell programs contained in references 1 and 3. The program in reference 1 is 
based on Sanders' shell theory for small displacements with moderate rotations (ref. 11). 
The governing equations of motion are reduced to six partial differential equations with 
first-order spatial derivatives and second-order time derivatives and are applicable to 
shells loaded axisymmetrically. The time derivatives are represented numerically by 
Houbolt’s backward -difference expression as outlined in references 1 and 4 and the spatial 
derivatives are represented by central differences. The equations are linearized using a 
Newton-Raphson technique and solved using Gaussian elimination. Modifications to the 
equations contained in reference 1, for the inclusion of rotary inertia and transverse 
shear effects, are presented in appendix A. 
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Sanders' shell theory for small rotations is also used in reference 3 which deals 
with the asymmetric behavior of shells. The equations are written as four governing 
partial differential equations with second -order derivatives in both space and time. The 
Houbolt backward-difference scheme and central differences are again used to represent 
the time and space derivatives, respectively, and Gaussian elimination is used to solve 
the equations. Since backward differences are used to represent time derivatives, the 
method is dependent on the previous time history of the shell. In general, this method 
allows one to take moderately large time steps without detrimentally affecting the quality 
of convergence. For structures in which high frequencies dominate the response, very 
small time increments must be taken. The convergence criterion used in this report was 
to use a time increment such that decreasing the time increment size by a factor of 2 
would not affect the results more than 1 percent after 10 time steps. 


Modal Superposition 

The method used for the modal superposition calculations is described in appen- 
dix B. Basically the approach is to take the stiffness and mass matrices and modal char- 
acteristics supplied by a previously executed program and compute generalized stiffness 
and mass matrices. These are then used to form a set of uncoupled equations which are 
solved in closed form for a given set of initial conditions. All supplied modes are used 
to preclude any arbitrary omission of modes which could have significant contributions. 

The modal characteristics are supplied either by the program of references 7 and 8 
or by a program based on the method described in appendix C. A critical feature of this 
approach is to obtain the modal data accurately. A discussion of this feature is contained 
in appendix C. 


ILLUSTRATIVE EXAMPLES 

The two shell configurations shown in figures 1 and 2 were considered. In these 
problems a satisfactory uniform grid spacing has been utilized unless otherwise indicated. 

Cylindrical Shell 

The first configuration is a cylindrical shell with freely supported boundary condi- 
tions and with a uniform initial velocity. This shell is similar to the one with cutouts 
used in reference 5. Since the structure has a plane of symmetry perpendicular to the 
shell axis and since the initial velocity is uniform, only half the structure needs to be 
analyzed. The boundary conditions at the initial edge (s = 0) are 
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N x = ° 

W = 0 
M x = 0 

and at the plane of symmetry (s = L/2) are 

U = 0 
Q = 0 
0 = 0 

The initial conditions were taken to be 

U(s,0) = 0 
V(s,0) = 0 
W(s,0) = 0 
0(s,O) = 0 

and 

U(s,0) = 0 
V(s,0) = 0 

W(s,0) = 2.54 cm/sec (1.0 in. /sec) 

/3(s,0) = 0 

where the displacements are all functions of space and time and a dot represents a deriv- 
ative with respect to time. The material properties were taken to be 

E = 6464 N/m 2 (o.9375 lb/in 2 ) 
p = 10.7 Gg/m^ (l.O lb-sec 2 /in 4 ) 
v = 0.25 
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The dimensions of the cylinders are characterized by the length-radius ratio L/r and 
the thickness -radius ratio h/r. Calculations are carried out for a range of L/r from 
0.3 to 2.4 and a range of h/r from 0.0125 to 0.10. 

The finite -element grid used in the modal superposition method is shown in fig- 
ure 3(a). For the direct integration method the spacial -difference increments used were 
As = L/96 and the nondimensional time -difference increment was taken to be 0.10 sec. 
The loading along the shell meridian midsurface is an inertial loading corresponding to 
the initial conditions on velocity. 


Conical Shell 

The second configuration was a cone similar to the Viking aeroshell (ref. 12). The 
finite -element grid of figure 3(b) is applied consistent with the geometry of figure 1(b). 
Here the dimensions are 

r A = 80.0 cm (31 .5 in.) 
rg = 164.0 cm (64.6 in.) 

L = 89.46 cm (35.22 in.) 

\l/ = tt/ 9 (20°) 
h = 0.101 cm (0.04 in.) 

and the material properties are 

E = 73.08 GN/m 2 (l .06 x 10 7 psi) 
v = 0.315 

p = 2713 kg/m 3 (2.54 x 10" 4 lb -sec 2 / in 4 ) 

The boundary conditions are freely supported at the inner edge (s = 0) and free at 
outer edge. Therefore, at s = 0: 

N x =0 

w = 0 

V = 0 
Mx ~ 0 
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and at s = L: 


N x = 0 

Q x = ° 

N x0 = 0 
M x = ° 

The initial conditions were 

U(s,0) = 0 
V(s,0) = 0 
W(s,0) — 0 
/3(s,0) = 0 

and 

U(s,0) = 0 
V(s,0) = 0 

W(s,0) = -v Q (v Q = 2.54 cm/see (1.0 in. /sec)) 

/3(s,0) = 0 

The finite -difference grid sizes were L/100 in the spacial direction and 10 -D sec in 
time increments. The nonhomogeneous, initial condition is a linearly varying normal 
velocity along the shell meridian with peak amplitude of v 0 at s = 0. This provides 
the only loading on the structure. 

RESULTS AND DISCUSSION 

These results are presented for future comparisons with large general-purpose 
shell codes. Brief parametric studies are presented as well as a tabular listing in 
appendix D of all the numerical results. 
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Cylindrical Shell 


Results for the cylindrical shell with loading corresponding to an initial velocity 
normal to the shell middle surface are shown in figure 4. . Both the direct integration 
and modal superposition approaches are compared with an "exact” solution similar to 
that described in appendix C of reference 13. Here, the exact solution is again a modal 
superposition technique where all the modes were obtained as sine and cosine waves. 

The series summation in the solution was truncated after 100 modes. The curves in fig- 
ure 4 represent the normal displacement W/r along the shell meridian at 1.2 seconds 
(after 12 steps in time for the direct integration scheme have been taken). The agree- 
ment is quite good except for the peak displacement region around s/L = 0.10 for the 
modal superposition approach; even then the discrepancy is only about 5 percent and this 
is attributed to the fact that only 10 elements were used to model the shell. The accuracy 
of the modal superposition approach is obviously dependent on obtaining a large number 
of reliable eigenvalues and eigenvectors. As the number of elements is increased, com- 
puter storage and time increase proportionally. Once the modal data are obtained, how- 
ever, the modal superposition approach can be applied to any point independent of time 
history; in fact, repeated solutions at various points in time can be obtained quite effi- 
ciently. Since the modal data are dependent only upon the shell geometry and boundary 
conditions, the modal superposition approach allows an additional degree of flexibility in 
efficiently studying a variety of initial conditions. 

Figures 5 to 8 show the effect of transverse shear deformation and rotary inertia 
(secondary effects) on transient response of the cylindrical shell for three values of L/r. 
The direct integration scheme is used to produce these results. In these figures the 
ratio h/r is held to a value of 0.05. In these plots the linear solution and the linear 
solution with the secondary effects included are compared. Figure 5 is a plot of normal 
displacement W/r as a function of the meridional distance at a time of 1.2 sec (12 steps 
in time). For this curve the secondary effects affect the solution only slightly; however, 
the moment resultant M x , which has a more prominent boundary region as shown in fig- 
ure 6, indicates that these terms have a more pronounced effect. 

Figures 7 and 8 show that these effects also grow as the shell is shortened while 
the other geometrical data remain constant. Since in these figures there is no pronounced 
boundary -region behavior, this deviation from the linear solution is attributed primarily 
to the transverse-shear-deformation terms. 

Figures 9 to 11 indicate the importance of being able to include nonlinear terms in 
the solution as well as the effect of varying the ratio h/r where the ratio L/r remains 
at 2.4. In these cases the direct integration scheme must be used since modal super- 
position is restricted to linear problems. Figure 9 is a comparison of the normal 
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displacement W/r along the meridian for the linear, nonlinear, and secondary effects 
solutions. The curves are shown at 2.4 sec (24 increments in time) for the cylinder. 

The secondary effects are of little importance as also shown in figure 5, but the non- 
linear terms greatly affect the character of the response shape. Figures 10 and 11 show 
that for thinner shells these terms remain just as large and tend to damp out the response 
in the region around the shell edge at s = 0. 

Figure 12 shows the results for the linear asymmetric response for the second 
harmonic (n = 2) when L/r is 2.4 and h/r is 0.05. In this figure the normal displace- 
ment is plotted as a function of meridional distance s/L at a time of 2.4 see (24 incre- 
ments in time for the direct integration scheme) for both approaches. Here the agree- 
ment of the methods is better than that in figure 4. In this case 20 elements along the 
meridian were taken. This better modeling improves the quality of the results; however, 
to obtain this accuracy all the modal data must be stored and as such may make this 
procedure impracticable in large-scale applications. 

In figure 13 both maximum normal displacement and center displacement are 
plotted as a function of time for the structure when L/r is 2.4 and h/r is 0.10. The 
solutions give a profile of time history for linear, nonlinear, and secondary effects for 
axisymmetric analyses. 

Conical Shell 

The conical shell shown in figure 1(b) is analyzed for a load corresponding to a 
linearly decreasing initial normal velocity along the shell meridian. The results in fig- 
ure 14 are for the normal displacement corresponding to n = 2 at a time of 
20 x 10" 6 sec (or after 20 time increments in the direct integration scheme). Satis- 
factory agreement is obtained by using a fine finite -element grid (e/L = 0.025) in the 
region near the smaller boundary of the shell and coarser grids with larger elements 
throughout the remaining portion of the shell (see fig. 3(b)). The direct integration 
approach gives an almost linear decay of the displacement to zero at the outer edge 
(s/L = 1.0). The small oscillatory character of the modal solution is attributed to the 
coarsest grid size (e/L = 0.100) used near the larger boundary of the shell. 

The two approaches each offer unique capability. The modal approach yields 
numerical accuracy at any point in time. Also, the modal data are obtained independent 
of the initial conditions. Thus, the modal superposition results may be obtained at any 
time point and for any admissible set of initial conditions. 

The distinctive features of the direct integration approach are that high-frequency 
response may be analyzed and that both secondary and nonlinear effects can be easily 
included. 
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CONCLUDING REMARKS 


Both direct integration and modal superposition schemes for transient response 
analysis of shells of revolution have been discussed. The analysis reported in NASA 
TN D-6158 has been extended to include the secondary effects of rotary inertia and trans- 
verse shear deformation terms. In addition, the analytical scheme for modal superposi- 
tion and a finite -element analysis for obtaining the axisymmetric mode shapes and fre- 
quencies of cylindrical shells with secondary effects included have been presented. 

A brief parametric study of nonlinear and secondary effects has been shown and 
tabular data of the results have been made available. In general, the inclusion of non- 
linear terms is of much more importance than the inclusion of secondary effects although 
there are thin-shell structures where these secondary effects will be important. The 
advantages of the direct integration approach are that geometric nonlinear and secondary 
effects are easy to include and high-frequency response may be calculated. In compari- 
son to the modal superposition technique the computer storage requirements are less. 

The advantages of the modal superposition approach are that the solution is independent 
of the previous time history and that once the modal data are obtained, the response for 
repeated cases may be efficiently computed. Also, any admissible set of initial conditions 
can be used with the modal superposition method. The actual preference of one method 
over the other is dependent, therefore, on the problem requirements. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Hampton, Va., January 28, 1974. 
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APPENDIX A 


MODIFICATIONS OF ANALYSIS OF NASA TN D-6158 TO INCLUDE 
ROTARY INERTIA AND TRANSVERSE SHEAR DEFORMATION 

To include rotary inertia and shear, equations (4) and (5) of reference 1 should be 
replaced by 

m 'u + £ 2 r i ( m n - m 22) -p + n^ = 0 <A1) 

2(1 + V i n) 

-1_ l £1 q = w * - p u - 0 (A2) 

k E 10 h 


where q is the transverse shear resultant and 0 the rotary inertia. These changes 
require that equations (A6) and (A9) in reference 1 be changed to include, respectively, 


h 52 


2 ( 1 + y 12 ) 

Ejollk 2 


(A3) 



The notation of equations (Al) to (A4) is the same as that used in reference 1 


(A4) 
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APPENDIX B 


MODAL SUPERPOSITION PROCEDURE 

This appendix describes the procedure based on reference 6 for obtaining the 
transient response of shell structures using modal superposition for problems with non- 
homogeneous initial conditions. In this analysis for any harmonic number n, the 

, the eigenvector 
. ., m are provided by an 
analysis such as that in reference 8 or that described in appendix C. The equations of 
motion from reference 6 may be written as 

[ M n]{y n ) + [ K n]{ y n] = 0 < B1 ) 


m x m stiffness matrix |K n J ? the m x m mass matrix 
matrix anc * ^he eigenvalues X* where i - 1, 2 , . 


where the dots indicate differentiation with respect to time. The initial conditions are 


w = w 

{*„} = {y 0 } 

Transforming to the modal coordinates |q n j g iveS 

W = NW 

where the ith column of |V , denoted as satisfies the equation 

[ k JK}=4K]R} 


(B2) 


(B3) 


(B4) 


with xj^ being the eigenvalue corresponding to the eigenvector The eigenvectors 

are normalized so that the generalized mass for each mode is unity 


M T K][*n]=H 

Using equation (B5) with equation (B4) gives the additional result that 

W T NK>['n] 


(B5) 


(B6) 



APPENDIX B 


where is a diagonal matrix of the eigenvalues. Using equations (B3), (B5), 

and (B6) in equation (Bl) gives a set of equations for the elements of jq n j 

{«n} + [%]{%) = ° 


The initial conditions (eqs. (B2)) expressed in terms of |q n J are 

Mo=M T MM 1 

Wo = M T MW j 


(B7) 


(B8) 


The solution of equation (B7) subject to the conditions of equations (B8) is easily obtained. 
The ith element of fqA denoted as is given by 



sin 


X* t + b* cos 
n n 



(B9) 


where i = 1, 2, . , ., m. 


Then 



(BIO) 

(Bll) 




(B12) 


The required result |y n J is finally obtained by substituting equation (B12) into 
equation (B3). 
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APPENDIX C 


FINITE -ELEMENT CALCULATION OF AXISYMMETRIC VIBRATION 
MODES OF A CYLINDER INCLUDING SECONDARY EFFECTS 


In figure 2 the displacement components in the meridional and normal directions at 
a general point in the cylinder cross section are denoted as U and W, respectively. 
The corresponding displacement components at a point on the middle surface of the 
cylinder are denoted by u and w, respectively. The rotation of a shell generator is 
denoted by 0. A point is located in the cylinder by the meridional coordinate s, the 
circumferential coordinate 0, and the normal coordinate z. The middle surface of the 
shell is taken as the reference surface and is located a distance r from the shell axis. 
Following the usual convention, the U and W, at any time t, for axisymmetric vibra- 
tion at frequency o> may be separated into functions of displacement and time as 
follows: 


U(s,z,t) = [u(s) + z/3(s)]e ia,t 
W(s ,z,t) = w(s)e lcot 


(Cl) 


The maximum strain energy is 


+ 02 e 2 + a^y^r ds d0 dz 


The strain -displacement relations are 



(C3) 
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APPENDIX C 


The constitutive equations are 


CT 1 = 
a 2 = 

a 13 


-^r( e 

1 - v 2 v 


= k 2 Ghy 13 


+ ke 2 ) 

+ 1 ' e l) 






(C4) 


Using equations (Cl), (C3), and (C4) in equation (C2) and performing the integration 
with respect to the variables 6 and z give the following form for the maximum strain 
energy: 


where 



2v * 

+ — u f w 
r 



+ D (/ 3-) 2 



(C5) 


and 



D = 


Eh 3 

12(l - v 2 ) 


The integration indicated by \ is taken over the cylinder length. The membrane stiff- 

J s 

ness is C and the flexural stiffness is D, Equation (C5) may be written as a quadratic 
form as follows: 


v = f J s {«} T [R]{a}d S 

where 


[R] = 


2Cr 

2i^C 

0 

0 

0 


2vC 

2C/r 

0 

0 

0 


0 

0 

2k 2 Ghr 

0 

0 


0 

0 


2k 2 Ghr 


0 

0 

0 

0 

2Dr 


(C6) 


(C7) 
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w 


{6} = < w’} 


0 


L^J 


(C8) 


Kinetic Energy 

The maximum kinetic energy of a thin cylindrical shell is 


T = | JJJ p[u 2 + W^jr ds d6 dz 


(C9) 


Substituting equation (Cl) into equation (C9) and carrying out the integration with 
respect to 6 and z yield the following result for the maximum kinetic energy: 


T = u>‘ 


>J 7r J phr(u 


2 + w 2 + a/3 2 )d& 


(CIO) 


where a = Equation (CIO) may be written as a quadratic form as follows: 


(Cll) 


where 


M - 


2phr 

0 

0 


0 

2phr 

0 


0 

0 

2aphr 


(C12) 
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Representation of Shell Finite Elements 


(Cl 3) 


The cylinder is modeled by finite elements as shown in figure 3. Each element is 
an exact slice of the cylinder being analyzed. The following definitions are made in con- 
nection with the finite -element representation: 
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K total number of elements 

e k length of the kth element, k - 1, 2, . . K 

x meridional coordinate inside kth element measured from center of element; 



s k distance from first edge of cylinder to center of kth element 

From these definitions for x and s k it follows that 
s = s k + x 

The convention used herein is that quantities such as displacement components and 

e k e k 

their derivatives at s = s k - an( j s + _ are denoted by the subscripts k 
and k + 1, respectively. For example, u k is the meridional displacement at 

^k t ^k 

s = s k - and 0 k+ j is the derivative of /3 at s = s k + — . 

The quantities u, w, and ft are approximated over each element by the following 
polynomials: 

u=a 0,k + a l,l^ + a 2,k x2 + a 3,l ! x3 
w = b 0,k + b l,k x+b 2,k x2+b 3,k x3 
^ = c 0,k + c l > k x + c 2,k x2 + c 3,k x3 

Using equations (C8) and (C14) yields 

( 5 >= [ X ]<V) (CIS) 
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where 



(Cl 6) 


M- 


0 1 2x 3x^ 


0 0 

X x^ 


0 0 


0 0 
„3 


1 2x 3x 


0 0 0 
0 0 0 0 


■, 2 3 

1 XXX 


0 1 2x 3x^ 


(C17) 


“ 6 , 


T . « C Ir C lr 

Inserting x = -y and y into appropriate locations in equation (C17), using 
equation (Cl 5) with the first of equations (Cl 4) , and inverting the resulting equation give 


M = K] W 


(CIS) 
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where 



(C19) 
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Formulation of Element Stiffness Matrix 

Substitution of equation (Cl 5) into equation (C6) and integration over a single ele- 
ment give 


‘H T NM 


(C20) 




(C21) 


The integral in equation (C21) is evaluated numerically by using the trapezoidal rule with 
100 integration intervals. Using equation (C18) in equation (C20) gives 


v = 5W T M T N[ A k]W 


(C22) 


Thus the element stiffness matrix for the kth element is given by 


K] - [ A k] T [ c k][ A k] 


(C23) 


Formulation of Element Mass Matrix 


Using equations (Cl 3) and (C14) gives 


where 


to-HM 


y] = o 


(C24) 


1 x x 2 x 3 0 0 000000 


0 0 0 1 x x 2 x 3 0 0 0 0 


[o 000 0 0 0 0 1 x x 2 x 3 

(C25) 

Substitution of equation (C24) into equation (Cll) and integration over a single ele- 
ment yield the following expression for maximum kinetic energy: 




(C26) 
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where 

The integral in equation (C27) is evaluated by use of the trapezoidal rule. Using equa- 
tions (C18) and (C26) gives 

Thus the element mass matrix is given by 

[ M kJ - [ A k] T [ z k][ A k] < c29 > 

Formulation of Modal Equations 

The stiffness and mass matrices for the complete cylinder are synthesized from 
the corresponding element matrices as described in reference 7. This synthesis is 
accomplished by superimposing the lower right 6x6 block of the kth matrix and the upper 
left 6x6 block of the k+lst matrix. The superposition of element matrices is predicted 
on the following conditions of compatibility at junctures between elements: 


dx 


(C27) 


r -\ 
u k+l 


r *\ 

u k+l 

w k+l 


w k+l 

^k+1 

Vu' 

u k+l 

, .« 

%+l 1 
u k+l 

w k+l 


w k+l 

Iv. 

kth element 

M: 


for all k < K 


(C30) 


k+lst element 


By denoting the master stiffness and mass matrices for the shell by |J3j and 
respectively , and using the principle of minimum potential energy one obtains 



> 



0) 


M 


w-° 


(C31) 


which is the modal set of equations for the system and where |yj is a 6(K+1) x 1 
vector given by 
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<?>= 



(C32) 


The matrices in equation (C31) are modified by deleting rows and columns as necessary 
to satisfy the specific edge constraints. Equation (C31) is solved by a standard eigen- 
value extraction algorithm to yield natural frequencies co and modal columns {y 


Mode Shapes and Modal Stress and Moment Resultants 

The modal columns {y} contain values of u, w, /3, u' , w’, and /3' at the 
element junctures. In order to obtain a detailed mode shape, the coefficients of the poly- 
nomial representations of u, w , and are obtained by use of equation (Cl 8). Then, 
by using equation (C14), the mode shapes may be computed within each element and hence 
over the entire shell. 

The modal stress and moment resultants are obtained from the modal strains and 
displacements. (See appendix B of ref. 8.) 

Meridional stress resultant: 

N x = C ( u ’ + t ) (C33) 

Meridional moment resultant: 

M x = D/3' (C34) 


Transverse shear stress resultant: 

Q x = k 2 Gh(w' + j3) (C35) 

Accuracy of Modal Characteristics 

To assess the reliability of the modes supplied by this method, a computer program 
based on the method was used to calculate mode shapes, frequencies, and modal stresses 
of a cylinder whose properties are listed below: 
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L = 63.5 cm (25 in.) 


r = 25.4 cm (10 in.) 
h= 1.01 cm (0.4 in.) 

E = 207 GN/m 2 (3 x 10 7 lb/in 2 ) 
p = 7843 kg/m 3 (7.33 x 10" 4 lb-sec 2 /in 4 ) 


v = 0.30 
k 2 = 0.87 


Calculations are carried out for clamped -clamped boundaries. Results were com- 
pared with those obtained by the more general-purpose computer program (ref. 7). The 
latter program does not include rotary inertia or transverse shear effects but a satis- 
factory comparison can be made for the low-frequency modes where effects of rotary 
inertia and transverse shear deformation are small. 

There is a certain amount of difficulty in deciding which modes to compare. For 
the low-frequency modes, close agreement is observed between the two methods. 

Because of the differences noted in the two programs, there is not a one-to-one corre- 
spondence between modes from each program over the entire frequency spectrum. The 
following table contains frequency comparisons for the 10 lowest frequencies. The fre- 
quencies are in good agreement for those modes, with an increasing disparity evident 
with increasing mode number. This disparity in the higher frequencies is due to the 
secondary effects. 


COMPARISON OF FREQUENCIES FROM PROGRAM METHOD OF APPENDIX C 


AND REFERENCE 7 FOR A CLAMPED -CLAMPED CYLINDER 


Description of 
mode 

Frequency, Hz 
(appendix C) 

Frequency, Hz 
(ref. 7) 

First w mode 

3120.38 

3122.17 

Second w mode 

3259.15 

3262.10 

Third w mode 

3367.92 

3370.17 

Fourth w mode 

3410.42 

3420.52 

Fifth w mode 

3700.37 

3722.45 

Sixth w mode 

4045.63 

4085.35 

First u mode 

4465.32 

4472.51 

Seventh w mode 

4644.45 

4725.28 

Eighth w mode 

5363.40 

5497.80 

Ninth w mode 

6203.20 

6414.83 


25 




APPENDIX C 


Mode shapes and modal stress comparisons for the first five modes are shown in 
the machine plots in figure 15. On the left side of each frame are results from the pres- 
ent program and on the right are results from the program of reference 7. Since rotary 
inertia and transverse shear deformation are neglected in reference 7, there are no cor- 
responding plots for the rotation 0 or the transverse shear resultant Q x in the right 
side of the figure. The corresponding curves in the two analyses for the displace- 
ments u and w have excellent agreement in mode shape and, further, the maximum 
values are in agreement through at least two and often three significant figures. Also 
modal stress and moment resultants are in excellent qualitative agreement. There is 
agreement of the first significant digit in peak stress and moment resultants between the 
two programs and in many cases there is agreement through three significant figures. 

The oscillations in the plot for Q x are at locations which correspond to element 
junctures and are attributed to the fact that continuity of the slope at element junctures 
of Q x is not enforced. Cusps are also evident in the meridional stress resultant in the 
fifth mode for the same reason. 
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TABULAR DATA 

Tables I to EX comprise listings of all data contained in the body of the report. 
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********************************************* ************************************ 

* TABLE I * 

* DATA PLOTTED IN FIGURE 4 FOR CYLINDER AT TIME *1.20 SEC * 

* N = 0 * 

* L/R *2.40 * 

* H * .0500 * 

************************************* ******************************************** 

* " * 


* 

X 

EXACT 


DIR 

I NT 

MOD 

SUP 

* 

* 

* 

S/L 

W DISP 

S/L 

W DISP 

S/L 

W DISP 

T 

* 

* 

* 

0*0000 

0.000000 

0.0000 

0.000000 

0.0000 

.000000 

* 

* 

.0100 

. 130890 

.0104 

. 135972 

. 0100 

.123230 

* 

* 

.0200 

.257700 

.0208 

.270025 

.0200 

.245474 

* 

* 

.0300 

.381580 

.0313 

. 400230 

.0300 

.364561 

* 

* 

.0400 

.502780 

.0417 

. 524646 

.0400 

.477949 

* 

* 

.0500 

.615877 

.0521 

.641326 

.0500 

.583886 

* 

* 

.0600 

.717844 

.0625 

.748341 

.0600 

.681938 

* 

* 

.0700 

.811291 

.0 7^9 

.843823 

.0700 

.770777 

* 

* 

. 0800 

.895611 

.0833 

.926019 

.0800 

.84 8469 

* 

* 

.0900 

.964090 

.0938 

.993359 

.0900 

.913440 

* 

* 

.1000 

1.015938 

.1042 

1.044555 

.1000 

.965138 

* 

* 

. 1100 

1.055888 

.1146 

1.078728 

.1100 

L. 004190 

* 

* 

. 1200 

1.082960 

.1250 

1.095597 

.1200 

1.030655 

* 

* 

. 1300 

1.091880 

.1354 

1.095717 

.1300 

1. 044593 

* 

* 

. 1400 

1.085650 

.1458 

1.080737 

. 1400 

1.046735 

* 

* 

. 1500 

1.070193 

.1563 

1.053573 

.1500 

1.038839 

* 

* 

. 1600 

1.046571 

.1667 

1.018386 

. 1600 

1.023836 

* 

* 

.1700 

1 .015252 

.1771 

.980244 

.1700 

1.004532 

* 

* 

. 1800 

.980116 

.1875 

• 94446 7 

. 1800 

.982871 

* 

* 

.1900 

.945075 

.1979 

.915738 

.1900 

.960563 

* 

* 

. 2000 

.915575 

.2083 

.897238 

.2000 

.939787 

* 

* 

.2100 

.896207 

.2188 

.890047 

.2100 

.923622 

* 

* 

.2200 

.884001 

.2292 

.893040 

.2200 

.913941 

* 

* 

.2300 

.878456 

.2 396 

.903316 

.2300 

.9 10017 

* 

* 

.2400 

. 883512 

.2500 

.917054 

.2400 

.909999 

* 

* 

.2500 

.899628 

.2604 

.930510 

.2500 

.912660 

* 

* 

.2600 

.919119 

.2708 

.940865 

.2600 

.918421 

* 

* 

.2700 

.938611 

.2813 

• 946 703 

.2700 

.926491 

* 

* 

.2800 

.951643 

.2917 

.948033 

.2800 

.934312 

* 

* 

.2900 

.958390 

.3021 

.945930 

.2900 

.939288 

* 

* 

.3000 

.960562. 

.3125 

.941982 

.3000 

.940423 

* 

* 

. 3 100 

.957999 

.3229 

.937737 

.3100 

.939391 

* 

* 

. 3200 

.947039 

.3333 

.934318 

.3200 

.938046 

* 

* 

. 3300 

.932425 

.3438 

.932259 

.3300 

.936989 

* 

* 

. 3400 

.918447 

.3542 

.931555 

.3400 

.935930 

* 

* 

. 3500 

.910884 

.3646 

.931854 

.3500 

.934033 

* 

* 

.3600 

.914963 

.3750 

.932668 

.3600 

.930924 

* 

* 

. 3700 

.924435 

.3854 

.933560 

.3700 

.9284 74 

* 

* 

. 3800' 

.934705 

.3958 

.934243 

.3800 

.927999 

* 

* 

.3900 

.941171 

.4063 

.934607 

.3900 

.929542 

* 

* 

. 4000 

.941697 

.4 16 7 

.934674 

.4000 

.931959 

* 

* 

.4100 

.939459 

.4271 

.934545 

.4100 

.933449 

* 

* 

.4200 

.936103 

.4375 

.934334 

.4200 

.933583 

* 

* 

.4 300 

.932394 

.4479 

.934132 

.4300 

.932939 

* 

* 

. 4400 

.928685 

.4583 

.933987 

.4400 

.932363 

* 

* 

.4500 

.926157 

.4688 

.933909 

.4500 

.932472 

* 

* 

.4600 

.92484 2 

.4792 

.933883 

.4600 

.933005 

* 

* 

.4700 

.927743 

.4896 

.933881 

.4700 

.932876 

* 

* 

.4800 

.932678 

.5000 

.933884 

.4800 

.931949 

* 

* 

.4900 

.937481 



.4900 

.930851 

* 

* 

.5000 

.940135 



.5000 

.930367 

* 


********************************************************************************* 


28 



************■****'* tt******* 


APPENDIX D 


************ ******* ****************** ** ***** ************ * ****************** * 

* TABLE II * 

* DATA FROM FIGURES 5, 6 FOR CYLINDER AT N = 0 ANO TIME = 1.2 SEC * 

* L/R = 2.40 * 

* H/R = 0.05 * 

** ********** *** ********** ** ********** *************************************** 


* * 

* W/R M-X * 

* * 

** **** ******************* ************ ** ** * ******************************* * * * 
* S/L LINEAR SEC EFF LINEAR SEC EFF * 

** **** ******* **** ******** ******** ******************* ** ********************** 


* 


o.oooo 

.0104 
. 0208 
.0313 
.0417 
.0521 
. 0625 
.0729 
. 0833 
.0938 
. 1042 
.1146 
.1250 
.1354 
. 1485 
.1563 
. 1667 
.1771 
.1875 
. 1979 
. 2083 
.2188 
.2292 
.2396 

* . 2500 

* .2604 

* .2708 

* .2813 

* .2917 

* .3021 

* .3125 

* .3229 


0.000000 

0.000000 

0. 

0. 

* 

.135972 

. 137643 

3.1934E-05 

3.2218E-05 

* 

.270025 

.273384 

6.4058E-05 

6 • 4735E-05 

* 

.400230 

.405260 

9.6456E-05 

9.7709E-05 

* 

.524646 

.531284 

L.2900E-04 

1 • 3097E-04 

* 

.641326 

.649437 

1.6130E-04 

1.6393E-04 

* 

.748341 

.757711 

1.9268E-04 

1.9577E-04 

* 

.843823 

. 854174 

2. 2220 E— 04 

2.2609E-04 

* 

.9260 L 9 

•937009 

2.4871E-04 

2.5523E-04 

* 

.993359 

1.004432 

2.7072E-04 

2.8311E-04 

* 

1.044555 

1.054609 

2.8613E-G4 

3.0647E-04 

* 

1.078728 

1.085779 

2.9187E-04 

3. 1685 E- 04 

* 

1.095597 

1.096842 

2. 83 76F-04 

3.0215E-04 

* 

1.095717 

1.088320 

2.5718E-04 

2. 5267E-04 

* 

1.080737 

1. 063213 

2.0858E-04 

1.6833E-04 

* 

1.053573 

1.027130 

1. 3 7 89 E— 04 

6. 2092E-05 

* 

1.018386 

. 987378 

5.0530E-05 

-4.3294E-05 

* 

.980244 

.95 1254 

-4.2020E-05 

-1.2389E-04 

* 

.944467 

.924248 

-1. 2406 E— 04 

-1.63B7E-04 

* 

.915738 

.908876 

-U 79786-04 

-1.6076E-G4 

* 

.897238 

.904555 

-1.9831E-04 

-1.2468E-04 

* 

.890047 

.908427 

-1.7757E-04 

-7.2789E-05 

* 

.893040 

.916664 

-1.2541E-04 

-2 . 2302E-05 

* 

.903316 

.925727 

-5.72 06E-05 

1.4875E-05 

* 

.917054 

.933172 

9.0056E-06 

3.4206E-05 

* 

.930510 

.937889 

5 • 8039E-05 

3. 7530E-05 

* 

.940865 

.939890 

8. 1640E-05 

3 . 0322E-05 

* 

.946703 

.939860 

7. 9806 E- 05 

1 . 8765E-05 

* 

.948033 

.938698 

5 • 9327E- 05 

7.6974E-06 

* 

.945930 

.937189 

3.0394E-05 

-2.3200E-07 

* 

.941982 

.935846 

2. 8831E-06 

-4.3886E-06 

* 

.937737 

. 93489? 

- 1 * 6435E-05 

-5.4513E-06 

* 


* .3333 

* .3438 

* .3542 

* .3646 

* .3750 

* .3854 

* . 3 95 8 

* .4063 

* .4167 

* .4271 

* .4375 

* .4479 

* .4583 

* .4688 

* .4792 

* .4896 

* .5000 


934318 

.934361 

-2.5025E-05 

-4.6256E-06 

* 

932259 

.934140 

-2.4115E-05 

-3.G569E— 06 

* 

931555 

.934105 

-1.7197E-05 

-1.5348E-06 

* 

931854 

.934148 

-8. 2866 E- 06 

-4.4230E-07 

* 

932668 

.934197 

-5.7000E-07 

•U 50006-07 

* 

933560 

.934221 

4. 2506 E- 06 

3.4930E-07 

* 

934243 

.934214 

5.9341E-06 

3. 12606-07 

* 

934607 

.934182 

5. 23 77E- 06 

1.8100E-07 

* 

934674 

•934139 

3.3174E-06 

4. 5500E-Q8 

* 

934545 

.934095 

1.2311 E— 06 

-5. 1400E— 08 

* 

934334 

.934055 

-3 • 37 20E-Q7 

- 1 . 023GE-07 

* 

934132 

.934024 

-1. 1321E-06 

-1. 1680E-07 

* 

933987 

.934000 

-1.2303E-06 

— 1. 1020E-07 

* 

933909 

.933984 

-8.8670E-07 

-9.5320E-07 

* 

933883 

.933973 

-3. 95306-07 

-8.0600E-08 

* 

933881 

.933967 

7.0000E-10 

-7.0700E-08 

* 

933884 

.933965 

l • 4970E-Q7 

-6.7200E-08 

* 


** ***** ************ ****** ********** ***************************************** 
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******************* ********************************************************* 

* TABLE III * 

DATA FROM FIGUKES 7, 0 FOR CYLINDER AT TIME = 1.2 SEC * 

N ■* 0 * 

H/R * 0.05 * 

4c ******** **************^*** ******************************* ******** ******** 

♦ 

L/R = 1.2 L/R - 0.3 * 

* 

******************************************************************** ******** 

* S/L LINEAR SEC EFF LINEAR SEC EFF * 

**************************************************************************** 

* * 


* 

0.0000 

0.000000 

0. 00 CO 00 

0.000000 

0.000000 

* 

* 

.0104 

.068162 

.069017 

.018187 

.018883 

* 

* 

.0208 

.136087 

• 1378C0 

.036356 

.037758 

* 

* 

.0313 

.203534 

.206104 

.054490 

.056603 

* 

* 

.0417 

.270266 

. 273689 

.072571 

.075406 

* 

* 

.0521 

.336042 

.340311 

. 090583 

.094150 

* 

* 

. 0625 

.400621 

.405722 

.108507 

. 112822 

* 

* 

.0729 

.463760 

.469673 

. 126326 

.131407 

* 

* 

.0833 

.525216 

.531908 

. 144022 

.149889 

* 

* 

.0938 

.5 84 746 

.592170 

.161578 

. 168253 

* 

* 

. 1042 

.642105 

.650202 

. 178977 

.186484 

* 

* 

.1146 

.697051 

. 705753 

.196200 

.204565 

* 

* 

.1250 

.749343 

.758580 

.213231 

.222481 

* 

* 

.1354 

.798744 

.808454 

.230052 

.240215 

* 

4c 

. 1458 

.845024 

.855157 

.246646 

.257750 

♦ 

* 

.1563 

.887956 

.898480 

.262996 

•275068 

* 

* 

.1667 

.927326 

.938207 

.279083 

.292153 

* 

* 

.1771 

.962932 

.974104 

.294892 

.308987 

* 

* 

.1875 

.994584 

1.005915 

.310405 

.325552 

* 

* 

. 1979 

1.022116 

1.033352 

.325605 

.341828 

♦ 

* 

.2083 

1.045385 

1.05 6118 

.3404 75 

.357798 

♦ 

* 

.2188 

1.064279 

1.073924 

.355000 

.373443 

★ 

* 

.2292 

1.078729 

1.086535 

. 369162 

.388744 

* 

* 

. 2 396 

1.088718 

1.093815 

.382945 

.403680 

* 

* 

.2500 

1.094287 

1. 095772 

.396333 

.418233 

* 

* 

.2604 

1.095556 

1.092597 

.409311 

.432384 

* 

4c 

.2 708 

1.092729 

1.0846 86 

.421863 

.446113 

* 

* 

.2813 

1.086103 

1.072643 

.433975 

.459400 

* 

* 

.2917 

1.076075 

1.05 7254 

.445631 

.472226 

* 

* 

.3021 

1.063139 

1.039447 

.456818 

.484572 

* 

4c 

.3125 

1.047879 

i. 020234 

.467522 

.496419 

♦ 

4c 

.3229 

1.030951 

1.000638 

.477728 

.507748 

* 

4c 

.3333 

1.013063 

.961629 

.487425 

.518542 

♦ 

4c 

.3438 

.994938 

.964058 

.496600 

.528782 

* 

4c 

.3542 

.977283 

. 948606 

.505242 

.538450 

♦ 

4c 

. 3646 

.960750 

.935756 

.513338 

.547531 

♦ 

* 

.3750 

.945893 

.925774 

.520879 

.556009 

* 

4c 

.3854 

.933142 

.918718 

.527854 

.563868 

* 

4c 

.3958 

.922772 

.914454 

.534255 

.571095 

* 

4c 

.4063 

.914890 

.912689 

.540073 

.577675 

♦ 

4c 

.4167 

.909434 

.913015 

.545299 

.583597 

♦ 

* 

.4271 

.906188 

.914946 

.549928 

.588850 

* 

4c 

.4375 

.904799 

.917959 

.553952 

.593423 

* 

4c 

.4479 

.904822 

.921537 

.557367 

.597308 

♦ 

♦ 

.4583 

.905763 

.925196 

.560167 

.600497 

* 

4c 

.4688 

.907125 

.928511 

.562349 

.602984 

* 

♦ 

.4792 

.908461 

.931135 

.563910 

.604764 

* 

4c 

.4 896 

.909415 

.932814 

.564847 

.605833 

♦ 

4c 

.5000 

.909759 

.933391 

.565159 

.606190 

* 


**41 4c** ********************************************************************** 
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** ** *** ******* ******* ******** *** ******* ***** ******************* ************* 

* TABLE IV * 

* DATA FROM FIGURE 9 FOR CYLINDER AT TIME = 2.4 SEC AND N = 0 * 

* L/R = 2.40 * * 

* H/R * C. 050 * 

** ********* ******** ********** ****** ** ** ** *** * ** ** * ** * ***** ** *** **** ******* ** 


4c 

* 

* 


EXACT 


NONLINEAR 


* 

SEC EFF * 
* 


******************* ********************************************************* 

* S/L W DISP S/L W DISP S/L W DISPL * 

**************** ************************* **** ******************************* 

* * 


* 

0.0000 

0.000000 

0.0000 

0.000000 

0.0000 

0.000000 

* 

* 

.0 104 

.052062 

.0104 

. 063265 

.0104 

.052276 

* 

* 

.02 08 

.104591 

.0208 

.135995 

. 0208 

.105078 

* 

* 

.0313 

.158016 

.0313 

.224 706 

.0313 

.158870 

* 

* 

.0417 

.212681 

.0417 

.331067 

.0417 

.214037 

* 

* 

.0521 

.268811 

.0521 

.450772 

.0521 

.270835 

* 

* 

.0625 

.326483 

.062 5 

.573600 

.0625 

.329357 

* 

* 

.0729 

.385594 

.0729 

.685894 

.0729 

.389502 

* 

* 

.0833 

.445845 

.0833 

.775264 

.0833 

.450954 

* 

* 

.0938 

.506730 

.0938 

.834951 

.0938 

.513169 

* 

* 

. 1042 

.567538 

.1042 

.865117 

.1042 

.575373 

* 

* 

.1146 

.627363 

.1146 

.870973 

. 1146 

.636575 

* 

* 

.1250 

.685130 

.1250 

.859836 

.1250 

.695592 

* 

* 

.1354 

.739632 

.1354 

.838793 

.1354 

.751088 

* 

* 

.1458 

.789576 

. 1458 

.813447 

. 14 58 

.801629 

* 

* 

.1563 

.833652 

.1563 

.787543 

.1563 

.845755 

* 

* 

.1667 

.870597 

. 1667 

. 763 175 

.1667 

.882065 

* 

* 

. 1771 

.899279 

.1771 

.741242 

. 1771 

.909311 

* 

* 

.18 75 

.918783 

. 18 75 

.721925 

.1875 

.926501 

* 

* 

.1979 

.928494 

.1979 

.705054 

.1979 

.933001 

* 

* 

.2083 

.928180 

.2083 

.690343 

.2083 

.928618 

* 

* 

.2188 

.918059 

.2188 

.677528 

.2188 

.913676 

* 

* 

.2292 

.898844 

.2292 

.666423 

.2292 

.889075 

* 

* 

.2396 

.871770 

.2396 

.656943 

.2396 

.856336 

* 

* 

.2 500 

.838566 

.2500 

.649093 

.25 00 

.817635 

* 

* 

.2 604 

.801404 

.2604 

.642938 

.2604 

.775796 

* 

* 

.2 708 

.762782 

.2 708 

.638567 

.2708 

.734201 

♦ 

* 

.2813 

. 725375 

.2813 

.636057 

.2813 

.696545 

* 

* 

.2917 

.691833 

.29 17 

.635436 

.2917 

•666418 

* 

* 

.3021 

.664555 

.3021 

.636661 

.3021 

.646740 

* 

* 

.3125 

.645452 

.3125 

.639582 

.3125 

.639159 

* 

* 

.3229 

.635726 

.3229 

.643873 

.3229 

.643595 

* 

* 

.3 333 

.635695 

.33 33 

.648915 

.3333 

.658113 

* 

* 

.34 38 

.644704 

.3438 

.653689 

.3438 

.679215 

* 

* 

.3 542 

.661149 

.3542 

.656834 

.3542 

.702530 

* 

* 

.3646 

.682620 

. 3646 

.657054 

.3646 

.723737 

* 

* 

.3 750 

.706181 

.3750 

.653852 

.3750 

.739463 

* 

* 

.3854 

.728730 

.38 54 

. 648231 

.3854 

.747901 

* 

* 

.3958 

.747416 

.3958 

.642756 

.3958 

.749009 

* 

* 

.4063 

.760023 

.4063 

.640558 

.4063 

.744248 

* 

* 

.4167 

. 765269 

.4167 

.643533 

.4167 

.735996 

* 

* 

.4271 

.762957 

.4271 

.650775 

.4271 

.726826 

* 

* 

.43 75 

.753966 

.43 75 

.658455 

.43 75 

.718863 

* 

* 

.4479 

.740072 

.4479 

. 661568 

.4479 

.713380 

* 

* 

.4583 

.723656 

.45 83 

.656695 

.45 83 

. 710683 

* 

* 

.4688 

.707336 

.4688 

.644109 

.4688 

.710261 

* 

* 

.4792 

.693599 

.4792 

.627941 

.4792 

.711122 

* 

* 

.4896 

.684473 

.4896 

.614442 

« 4 8 96 

.712185 

* 

* 

.5 000 

.681278 

.5000 

.609195 

.5000 

.712639 

* 


********************************************* ******************************* 
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**************************************************************************** 

* TABLE V * 

* DATA PROM FIGURES 10, 11 FOR CYLINDER AT TIME = 2.4 SEC * 

* N - 0 * 

* L/R - 2.40 * 

******************************************************* * ** ****************** 


H/R = 0.025 


H/R - 0.0125 


**************************************************************************** 

* S/L LINEAR NONLINEAR LINEAR NONLINEAR * 

**************************************************************************** 


* 

0.0000 

0.000000 

0. 000000 

0.0000 GO 

0.000000 

* 

* 

.0104 

.074326 

.112469 

.106086 

.174000 

* 

* 

. 0203 

. 149920 

.250327 

.215448 

.402987 

♦ 

* 

.0313 

.227821 

.419098 

.330139 

.646568 

* 

* 

.0417 

. 308629 

.597629 

.449960 

.808619 

* 

* 

.0521 

.392323 

.747703 

.571801 

.851468 

♦ 

* 

.0625 

.478139 

.838337 

.689544 

.821327 

* 

* 

.0729 

.564507 

.865831 

.794634 

.776025 

* 

* 

.0833 

.649073 

.849797 

.877391 

.739095 

* 

* 

.0938 

.728817 

.814626 

.928953 

.711846 

* 

* 

. 1042 

.800265 

• 776969 

.943590 

.691380 

* 

* 

.1140 

.859801 

.744190 

.920900 

.675737 

* 

* 

.1250 

.904060 

.717797 

.86 72 31 

.664065 

* 

* 

.1354 

.930372 

.696898 

.795642 

.656204 

* 

* 

. 1458 

.937220 

.680151 

.723852 

.652114 

* 

* 

.1563 

.924637 

.666500 

.670132 

.651206 

♦ 

♦ 

.1667 

.894495 

.655340 

.647881 

.652209 

* 

* 

. 1771 

.850592 

.646466 

.660553 

.653916 

* 

♦ 

. 1875 

.7984 74 

.639978 

.699343 

.655989 

* 

* 

.1979 

.744951 

.636134 

.745655 

.658695 

* 

* 

.2083 

.697281 

.635137 

.778600 

.661843 

* 

* 

.2188 

.662079 

.636876 

.784734 

.664375 

* 

* 

. 2292 

.644094 

.640696 

.764855 

.665470 

* 

* 

.2396 

.645082 

.645440 

.733195 

.665995 

* 

* 

.2500 

.663126 

.649840 

.708589 

.667506 

* 

* 

.2604 

.692718 

.653181 

.702876 

.668836 

* 

* 

.2708 

.725827 

.655739 

.714384 

.666747 

♦ 

* 

.2813 

.753843 

.658336 

.730994 

.663508 

* 

* 

.2917 

. 769943 

.660906 

.740028 

.667724 

* 

* 

.3021 

.771100 

.661420 

.737050 

.678770 

* 

* 

. 3125 

. 758946 

.657308 

.727117 

.679314 

* 

♦ 

. 3229 

.739043 

.649331 

.718763 

.660350 

* 

* 

. 3333 

.718749 

.643782 

.716782 

.644581 

♦ 

* 

.3438 

.704483 

.648524 

.719748 

.659247 

* 

* 

.3542 

.699534 

.664637 

.723052 

.691180 

* 

* 

. 3646 

.703271 

.682003 

.723529 

.700960 

* 


.3750 

.712022 

.685437 

. 721424 

.678105 

* 

* 

.3854 

.721055 

.667779 

.718945 

.646915 

* 

* 

.3958 

.726669 

. 638061 

.717740 

.628493 

* 

* 

.4063 

.72746 5 

.616437 

.717784 

.629406 

* 

* 

.4167 

.724359 

.619008 

.718097 

.648845 

* 

* 

.4 271 

.719605 

.645107 

.717957 

.677308 

* 

* 

.4375 

.715480 

.678224 

.717389 

.700847 

* 

* 

.4479 

.713327 

.698329 

.716805 

.707927 

* 

* 

.4583 

. 713294 

.694345 

.716476 

.694117 

* 

* 

.4688 

.714670 

.668813 

.716365 

.666786 

* 

* 

. 4792 

.716467 

.633988 

.716316 

.636958 

* 

* 

.4896 

.717861 

.604850 

.716258 

.611738 

♦ 

* 

. 5000 

. 718375 

.593539 

.716228 

.601029 

* 


*******************.********************************************************* 
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**#****##**##*#***#**###*#******###**#****#****#**#**#*#*#*********##******* 

* TABLE VI * 

* DATA FROM FIGURE 12 FOR CYLINDER AT N = 2 AND TIME = 2*4 SEC * 

* L/R = 2.40 * 

* H/R ^ 0.05 * 

******************* ****************** **************4t************************ 

* * 


* MODAL DIRECT 

* 


* 

* 


**************************************************************************** 

* S/L LINEAR S/L LINEAR * 

******************************************************************* ********* 

* * 


* 

0. 0000 

.000000 

0.0000 

-.000000 

* 

* 

.0100 

.110069 

.0052 

.057667 

* 

* 

.0200 

.220853 

.0156 

.173063 

4c 

* 

.0300 

.332183 

.0260 

.288625 

* 

* 

.0400 

.443397 

.036 5 

.404408 

* 

* 

.0500 

.55396 1 

.0469 

.520385 

* 

* 

. 0600 

.663900 

.0573 

.636420 

4c 

* 

.0700 

.773125 

.0677 

.752239 

* 

* 

.0800 

.881105 

.0781 

.867417 

4c 

* 

.0900 

.987186 

.0885 

.981360 

4c 

4c 

. 1000 

1.090911 

.0990 

1.093302 

* 

* 

. 1100 

1.191703 

.1094 

1.202311 

* 

$ 

. 1200 

1.287683 

.1198 

1.307305 

♦ 

* 

.1300 

1.378048 

.1302 

1.407073 

* 

* 

. 1400 

1.462927 

.1406 

1.500323 

♦ 

* 

. 1500 

1.542430 

.1510 

1.585723 

4c 

4c 

. L600 

1.615423 

.1615 

1.661973 

* 

4c 

.1700 

1.679213 

,1719 

1*727869 

4c 

❖ 

. 1800 

1. 732228 

.1823 

1.782390 

4c 

4c 

. 1 900 

1.774241 

.192 7 

1.824779 

4c 

4c 

.2000 

1.805987 

.2031 

1.854628 

4c 

4c 

.2100 

1.828787 

,2135 

1.871954 

4c 

* 

.2200 

1.843666 

.2240 

1.877264 

4c 

4c 

.2300 

1.850644 

.2344 

1.871593 

4c 

4c 

.2400 

1.849569 

.2448 

1 . 856524 

4c 

4c 

.25 00 

1.841163 

.2552 

1.834172 

* 

4c 

.2600 

1 .827602 

.2656 

1.807121 

4c 

* 

.2700 

1.811071 

.2760 

1.778317 

4c 

4c 

.2800 

1.79365] 

.2865 

1.750893 

4c 

4c 

.2900 

1.777222 

.2969 

1.727918 

4c 

4c 

.3000 

1.763097 

.3073 

1.712079 

4c 

❖ 

.3100 

1.751581 

.3177 

1.705328 

4c 

4c 

.3200 

1.742742 

.3281 

1.708547 

4c 

4c 

.3300 

1.738163 

.3385 

1.721327 

4c 

4c 

.3400 

1.739855 

.3490 

1.741925 

♦ 

4c 

. 3500 

1.748663 

.3 594 

1.767480 

4c 

4c 

. 3600 

1.762783 

.3698 

1.794459 

4c 

4c 

.3700 

1.779201 

.3802 

1.819278 

4c 

4c 

.3800 

1.796722 

.3 906 

1.838952 

4c 

4c 

.3900 

1.815139 

.4010 

1.851614 

4c 

4c 

.4000 

1.833819 

.4115 

1.856775 

4c 

4c 

.4100 

1.850103 

.4219 

1.855266 

4c 

4c 

.4200 

1.860021 

.4323 

1.848873 

4c 

4c 

.4300 

1.862291 

.4427 

1.839814 

4c 

4c 

.4400 

1.858081 

.4531 

1.830194 

4c 

4c 

.4500 

1.849928 

.4635 

1.821605 

4c 

4c 

.4600 

1.840083 

.4740 

1.814963 

4c 

4c 

.4700 

1.828884 

.4844 

1.810589 

4c 

4c 

.4800 

1.817606 

.4948 

1.808454 

4c 

4c 

.4900 

1.808673 

.5000 

1.808454 

4c 

4c 

.5000 

1.805036 



4c 


**************************************************************************** 
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**************************************************************************** 


* TABLE VII * 

* DATA FROM FIGURE 13A FOR CYLINDER AT N= 0 * 

* L/R = 2.4 * 

* H/R = 0.10 * 

##**##*####*#####******#*#**$##*********#********# **A******* **************** 

* * 

* LINEAR SEC EFF NONLINEAR * 

* * 


it*#****#*#***#** ##*##******#******* ***************** ****** ****************** 

♦ TIME W-MAX TIME W-MAX TIME W-MAX * 

************************* ** ** ** * *** $* ******* ******** * *♦* **** ***** **** ** *** ** 


♦ 

* 

0.0000 

0.000000 

0.0000 

0.000000 

0.0000 

0.000000 

* 

* 

* 

. 1 000 

.107519 

.1000 

.100648 

. 1000 

.104346 

* 

* 

.2000 

.215949 

.2000 

.203263 

. 2000 

.203334 

* 

* 

.3000 

.328021 

.3000 

.309262 

.3000 

.300908 

* 

* 

• 4000 

.438330 

.4000 

.416733 

.4000 

.394365 

* 

* 

.5 000 

.544708 

.5000 

. 522937 

.5000 

.483585 

* 

* 

.6 000 

.647014 

.6000 

.625444 

.6000 

.567785 

* 

* 

.7000 

.742572 

.7000 

.722552 

.7000 

.646371 

♦ 

* 

.8000 

.830545 

.8000 

.813057 

. 8000 

.718636 

* 

* 

.9 000 

.910239 

.9000 

.895943 

.9000 

.783986 

* 

* 

1.0000 

.981036 

1 . 0000 

.970272 

1.0000 

.842270 

* 

* 

1.1000 

1.042389 

1.1000 

1.035228 

l. 1000 

.891846 

♦ 

* 

1.2000 

1.093823 

1.2000 

1.090157 

1.2000 

.931921 

* 

* 

1.3000 

1.135194 

1.3000 

1.134572 

1.3000 

.962023 

* 

* 

1.4000 

1.166994 

1.4000 

1.168132 

1.4000 

.982065 

* 

* 

1.5000 

1.187990 

1.5000 

1.190622 

1. 5000 

.991570 

* 

* 

1.6000 

1.197978 

1.6000 

1.201952 

1.6000 

.990527 

♦ 

* 

1. 7 000 

1.196964 

1.7000 

1.202158 

1.7000 

.978762 

♦ 

* 

1.8000 

1.185464 

1.8000 

1. 191402 

1.8000 

.963455 

* 

* 

1.9000 

1.165167 

1.9000 

1.171299 

1.9000 

.956173 

♦ 

* 

2.0000 

I * 135928 

2.0000 

1.140870 

2.0000 

.943894 

* 

* 

2.1000 

1.096791 

2.1000 

1.100602 

2. 1000 

.926687 

♦ 

* 

2.2000 

1.045421 

2.2000 

1.051130 

2.2000 

.903616 

* 

* 

2.3000 

.982919 

2.3000 

.993176 

2.3000 

.876775 

♦ 

* 

2.4000 

.908478 

2.4000 

.927812 

2.4000 

.848017 

* 
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** **** ********** *** ************* ************ ****** ********** #** ** ** * * ** a$c ^ * ** 

* TABLE VIII * 

* DATA FROM FIGURE 13B FOR CYLINDER AT N= 0 * 

* L/R = 2.4 * 

* H/R = 0.10 * 

************************************************************ **************** 

* , * 

* LINEAR SEC EFF NONLINEAR * 

* * 

******* * **************** ********************* *************** **************** 
* TIME W-L/2 TIME W-L/2 TIME W-L/2 * 


***************************** **************** ******************************* 
* * 


* 

O.uOOQ 

0.000000 

0.0000 

0.000000 

0. 0000 

0.000000 

* 

* 

.1000 

.099834 

. 1000 

.099834 

. 1000 

.099834 

* 

* 

.2 000 

.198674 

.2000 

.198674 

.2000 

.198674 

* 

* 

.3000 

. 2 95540 

.3000 

.295540 

. 3000 

.295540 

* 

* 

.4000 

.389471 

.4000 

.389471 

.4000 

.389471 

* 

* 

.3000 

.479540 

. 5000 

.479537 

. 5000 

.479538 

* 

* 

.6 000 

.564868 

.6000 

.564847 

.6000 

.564844 

* 

* 

.7 000 

• 6 4462 9 

... 7000 

.644556 

.7000 

.644539 

* 

* 

.8000 

. 717986 

. 8000 

.717875 

.8000 

.717829 

* 

* 

.9 000 

.783986 

.9000 

. 7840 82 

.9000 

. 783986 

* 

* 

1 . 0000 

.841564 

1.0000 

.842523 

1.0000 

.842270 

* 

* 

1.1000 

.889945 

1.1000 

.892629 

1. 1000 

.891846 

* 

* 

1.2000 

.929353 

1.2000 

.933923 

1.2000 

.931921 

* 

* 

1.3000 

.961424 

1.3000 

.966042 

1.3000 

.962023 

* 

* 

1.4000 

.988579 

1.4000 

.988763 

1.4000 

.981972 

* 

* 

1.5000 

1.0 12323 

1.5000 

1.002056 

1.5000 

.991569 

* 

* 

1.6000 

1.031492 

1.6000 

1.006125 

1.6000 

.990527 

* 

* 

1.7000 

1.041782 

1.7000 

1.00142 1 

1.7000 

.978762 

* 

* 

1.8 000 

1.037093 

1.8000 

.988570 

1.8000 

.956608 

* 

* 

1.9000 

1.011992 

1.9000 

.968145 

1.9000 

.924782 

* 

* 

2.0000 

.963914 

2.0000 

.940311 

2.0000 

.884215 

* 

* 

2.1000 

.894069 

2.1000 

.904440 

2. 1000 

.835971 

* 

* 

2.2000 

.806857 

2.2000 

.858910 

2.2000 

.781362 

* 

* 

2.3000 

.708391 

2.3000 

.801310 

2.3000 

.722376 

* 

* 

2.4000 

.604917 

2.4000 

.729167 

2.4000 

.662495 

* 


**************************************************************************** 
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APPENDIX D 


*********************************** *************** ****************** ******** 

* TABLE IX * 

* DATA FROM FIGURE 14 FOR CONICAL SHELL * 

* TIME = 20 *£-06 SEC * 

* N =2 * 

**************************************************************************** 


MODAL 


DIRECT 


**************************************************************************** 
* S/L LINEAR S/L LINEAR * 

************** ****************** ******************************** ************ 


* 

0. 0003 

-1. 2620796-19 

0.0000 

3.1832 32E— 19 

* 

* 

. 0050 

-3. 657 528 E -06 

.0051 

-1. 4462 02 E-05 

* 

* 

.0100 

-7.344286E-06 

.0152 

-1. 9778236-05 

* 

* 

.0150 

-1. 16 3 078 E -05 

.0253 

- L. 873 338E-05 

* 

* 

.0 200 

- 1 • 546 3 3 OE- 05 

.0354 

- 1. 850952E— 05 

* 

* 

.0250 

- 1.824325E-05 

.045 5 

-1.8401 74E-05 

♦ 

* 

.0300 

- L. 92476 IE- 05 

.0556 

-1. 826713E-05 

* 

* 

.0350 

- 1. 9Q32Q5E-0 5 

.065 7 

- 1 • 812024E-05 

* 

* 

.0400 

- 1 . 848 854E- 05 

.0758 

-l. 7965 75E-05 

* 

* 

.0450 

- 1. 824 1 99E- 05 

.0859 

-1. 780554E-05 

* 

* 

.0 500 

-1.8343815-05 

.0960 

-1. 7640 8 IF- 05 

* 

* 

. 0550 

-1.8306 80 E - 05 

.1061 

- 1. 747 233 E-05 

* 

* 

.0600 

-1.819370 £-05 

,1 162 

-1. 73006 7^-0 5 

* 

* 

. 065 0 

-1. tt l l 59 7E-05 

.1263 

-1.7126 256-05 

* 

* 

.0700 

- 1 . 8060 16 £-05 

.1364 

- 1 .6949386-05 

* 

* 

.0750 

- 1. 7982966-05 

.1465 

— 1.6770 33E— 05 

* 

* 

.0800 

-1. 789206E-05 

.1566 

-1.65 89 29 E-05 

* 

* 

.0650 

-1. 7812366-05 

. 166 7 

- 1 • 64064 7E— 05 

* 

* 

. 0900 

-1. 7 74 0 18E-05 

.1768 

- 1. 6222 04E-05 

* 

* 

.0950 

-1. 766L05E-05 

.1869 

— I • 603 6 14E— 05 

* 

* 

. 1000 

-1. 75 71 166-05 

.1970 

-1.584 893 E-05 

* 

* 

.1050 

-1. 7485816-05 

.2071 

-1. 566052E- 05 

* 

* 

.1100 

-1. 740326E-05 

.2172 

- 1 . 547 1 03E-05 

* 

* 

.1150 

-1. 732003E-05 

.2273 

- 1 • 52805 76- 05 

* 

* 

. 1200 

- 1 . 723 5826- 05 

.2374 

-1.508 92 3 E - 05 

* 

* 

. 1 250 

-1. 7148846-05 

.247 5 

- 1.4 897 09E-05 

* 

* 

.1300 

-1. 705 7 02E-0 5 

.2576 

- 1.47042 IE-05 

* 

* 

.1350 

- 1 . 6 9 6 8 95 E - 0 5 

.2677 

— 1.45 i 068E-05 

* 

* 

.1400 

- 1 • 6884 97E-05 

.2778 

- 1 .431654E— 05 

* 

* 

. 1450 

-L. 6798916-05 

.2879 

- 1.412 1866-05 

* 

* 

. 1500 

-1.6706 76E - 05 

.2980 

- 1. 392667E— 05 

* 

£ 

. 1550 

-1.66L 3346-05 

.3081 

-1.373 1026-05 

* 

* 

. 1600 

- 1. 652 383E- 05 

.3 182 

- 1. 3534946-05 

* 

* 

.1650 

- 1 . 6 43 6 85 £-05 

.3283 

- 1 • 3338486—05 

* 

* 

. 1700 

- 1. 634807 E— 05 

.3 384 

-1.3141676—05 

* 

* 

. 1 750 

- L . 625449E-05 

.3485 

-1.2944526-05 

* 

* 

. 1800 

-1.61 5 899E-05 

. 3586 

-1.2747076-05 

* 

* 

. 1850 

-1.606 7446-05 

.368 7 

— 1. 2549346-05 

* 

* 

.1900 

-1.5978726-05 

. 3 78 8 

— 1. 235 13 5E- 05 

* 

* 

. 1950 

- 1. 588804E-05 

.3889 

- 1.21531 36-05 

* 

* 

.2000 

-1. 579261E-05 

.3990 

— 1.1954686—05 

* 

* 

.2050 

- 1. 5696 11E-05 

.4091 

- 1 • 1756 02E-05 

* 


.2100 

- 1 • 560 281 E-05 

.4192 

- 1. 155 7 176-05 

* 

* 

.2150 

- 1. 551 1 96E- 05 

.4293 

-1. 13581 5E-05 

* 

* 

.2200 

-I. 542028E-05 

.4394 

-1.1158966-05 

* 

* 

.2250 

-1.5325076-05 

. 4495 

-1.09596 IE-0 5 

* 

* 

.2300 

- 1 .522 7.3.3 E-05 

.4596 

-I. 0760136-05 

* 

* 

.2350 

-1. 5129436-05 

.469 7 

-1.0560506-05 

* 

* 

.2400 

-1.5031 2 6 E-05 

.4798 

-1.0360756-05 

* 

* 

.2450 

- 1. 4936 77E-06 

.4899 

-1.0160896-05 

* 

* 

.2 500 

- 1 * 486 2 08 E-05 

.5000 

-9.9609 096—0 6 

* 


************** ***************** ****** *************************************** 
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APPENDIX D 


*********************************** *************** ************************** 

* TABLE IX CON CL * 

* DATA FROM FIGURE 14 FOR CONICAL SHELL * 

* TIME = 20 *£—06 SEC * 

* N * 2 * 

******************* **************** ***************************************** 


* * 

* MODAL DIRECT * 

* * 

************************************* *************************************** 
* S/L LINEAR S/L LINEAR * 

************** ************************************************************** 


* 


* 


£ 

.2600 

- i* 477 102E-05 

.5101 

-9. 760828E-06 

* 

* 

. 2700 

- 1.4541 88 E^ 05 

.5202 

-9.56Q650E-06 

* 

* 

. 2800 

-1.4202 30E -0 5 

.530 3 

-9.360381E-06 

* 

* 

.2900 

-L.395L22E-05 

.5404 

-9. 160027E-06 

* 

* 

. 3000 

-1.3889756-05 

.5505 

— 8.95 9592E-06 

* 

* 

. 3100 

- L * 3 82 045E- 0 5 

.560 6 

-8. 75908 IE-06 

* 

* 

. 3200 

- 1. 357086E— 05 

.5707 

-8.558499E-06 

* 

* 

. 3300 

- 1. 323 129E-0 5 

.5808 

-8.357 850E-06 

* 

* 

.3400 

-1. 298519E-05 

.5909 

-8. 15 713 7E — 06 

* 

* 

. 3500 

-1. 291485E-05 

.6010 

-7.956364E-06 

* 

* 

.3600 

-1.284515E-G5 

.6111 

- 7. 755 5 3 5E — 06 

* 

* 

. 3700 

-i. 259658E-05 

.6212 

- 7. 554652 E-06 

* 

* 

. 3800 

- 1 . 225 389E-05 

.6313 

-7.3537 1 9E- 06 

* 

* 

.3900 

-l. 200520E-05 

.6414 

-7. 152 7 37E-06 

* 

* 

.4 000 

-i. 193342E-05 

.651 5 

-6. 95171 IE-06 

* 

* 

.4100 

-1. 18641 2t -0 5 

.6616 

- 6. 75064 IE-06 

* 

* 

. 4200 

-1. 1614035-05 

.6717 

-6.549531E-06 

* 

* 

.4300 

-1. 127076E-05 

.6818 

-6.3483 82E — 06 

* 

* 

.4400 

-1. 1020916-05 

.6919 

—6* 1471 96E-06 

* 

* 

.4 500 

- 1.094955E-Q5 

. 7020 

-5. 945976E— 06 

* 

* 

. 4600 

- 1.087 880E-Q5 

.7121 

-5. 744 72 2 E-06 

* 

* 

.4700 

- 1 . 062 8 13E- 05 

.7222 

-5.543437E-06 

* 

£ 

. 4800 

-1.028198E-05 

.7323 

-5 • 342 1 22E-06 

* 

* 

.4900 

- 1. 0029 70E-0 5 

.7424 

-5. 140 7 79E— 06 

* 

* 

. 5000 

-9. 9ol 7 97E-06 

. 752 5 

-4. 9394G8E-06 

* 

* 

.5200 

-9. 842923E-06 

. 762 6 

-4. 7380UE-06 

* 

* 

- 5400 

-9. 333085E-06 

. 7727 

-4.536590E-06 

* 


.5600 

-8. 60 79815-0o 

.7828 

-4. 335146E-06 

* 

* 

. 5800 

-8.0911 30E-06 

.7929 

-4 . 133679E-C6 

* 

* 

... 6 000 

-7.9773796-06 

.8030 

-3. 932 1 9 OE— 06 

* 

* 

.6200 

- 7. 860464E-06 

.8 131 

-3. 7306 8 1 E-06 

* 

* 

. 6400 

-7. 34 3 7 8 7E-06 

.8232 

-3.52915 3E- 06 

* 

* 

.6600 

-6. 619221E-06 

.8333 

- 3. 32 7 6 06 E-06 

* 

* 

.6800 

-6. 1 03 482 E- 06 

.8434 

-3.126 041 E-06 

* 

* 

. 7 000 

-5. 98599LE-06 

.8535 

-2 • 92 446 OE— 06 

* 

❖ 

. 7200 

-5.869312E-06 

.8636 

-2.722862E— 06 

* 

* 

. 7400 

-5. 352 8 1 6E- 06 

.8737 

-2.52 1 248E-06 

* 

* 

. 7600 

-4.626548E-06 

.3838 

-2.3196 19E-06 

* 

* 

. 7800 

-4. 1098906-06 

.8939 

-2.il/976E-06 

* 

* 

. 8000 

-3.992979E-06 

.9040 

-L.916320E-06 

* 

* 

. 8200 

-3. 875073F— 06 

.9141 

— 1 * 71465 IE — 06 

* 

* 

.8400 

-3.358155E-06 

.9242 

- 1. 51 2969E-06 

* 

* 

. 8600 

-2. 632091E-06 

.9343 

-1.31 1275E— 06 

* 

* 

. 8800 

-2. 1 1 4320E- 06 

.9444 

- 1 • 1095 70E-06 

* 

* 

.9000 

-1.996420E-06 

.9545 

- 9 . 0785 37E-07 

* 

* 

.9200 

-1. 830537E-06 

.9646 

-7.061278E-07 

* 

* 

.9400 

- 1. 361 560E-06 

.9747 

— 5. 0439 OOE-O 7 

* 

* 

. 9600 

-o. 3426206-07 

.9848 

-3.026 22 IE-07 

* 

* 

.9800 

-1. 196 1 91 E- 0 7 

.9949 

-1.009266E-07 

* 

* 

i • 0000 

1. 167211E-08 

1.0 000 

-3 • 9925 16E- 10 

* 


************** *********************** *************************************** 


37 



REFERENCES 


1. Stephens, Wendell B.: Computer Program for Static and Dynamic Axisymmetric Non- 

linear Response of Symmetrically Loaded Ortho tropic Shells of Revolution. NASA 
TND-6158, 1970. 

2. Stephens, Wendell B.; and Fulton, Robert E.: Axisymmetric Static and Dynamic 

Buckling of Spherical Caps Due to Centrally Distributed Pressures. AIAA J., 
vol. 7, no. 11, Nov. 1969, pp. 2120-2126. 

3. Stephens, Wendell B.; and Robinson, Martha P.: Computer Program for Finite- 

Difference Solutions of Shells of Revolution Under Asymmetric Dynamic Loading. 
NASA TN D-6059, 1971. 

4. Johnson, Donald E.; and Greif, Robert: Dynamic Response of a Cylindrical Shell: 

Two Numerical Methods. AIAA J., vol. 4, no. 3, Mar. 1966, pp. 486-494. 

5. Geers, Thomas L.; and Sobel, Lawrence H.: Analysis of Transient, Linear Wave 

Propagation in Shells by the Finite Difference Method. NASA CR-1885, 1971. 

6. Raney, J. P.; and Howlett, J. T.: A Modal Solution for Wave Propagation in Finite 

Shells of Revolution. J. Spacecraft & Rockets, vol. 8, no, 6, June 1971, pp. 650-656. 

7. Adelman, Howard M.; Catherines, Donnell S.; and Walton, William C., Jr.: A Method 

for Computation of Vibration Modes and Frequencies of Orthotropic Thin Shells of 
Revolution Having General Meridional Curvature. NASA TN D-4972, 1969. 

8. Adelman, Howard M.; Catherines, Donnell S.; Steeves, Earl E.; and Walton, 

William C., Jr.: User's Manual for a Digital Computer Program for Computing 
the Vibration Characteristics of Ring-Stiffened Orthotropic Shells of Revolution. 
NASA TMX- 2138, 1970. 

9. MacNeal, Richard H., ed.: The NASTRAN Theoretical Manual (Level 15). NASA 

SP- 22 1(01), 1972. 

10. Almroth, B. O.; Brogan, F. A.; and Marlowe, M. B.: Collapse Analysis for Shells of 

General Shape. Volume 1 - Analysis. AFFDL-TR-71-8, U.S. Air Force, Aug. 
1972. 

11. Sanders, J. Lyell, Jr.: Nonlinear Theories for Thin Shells. Quart. Appl. Math., 

vol. XXI, no. 1, Apr. 1963, pp. 21-36. 

12. Heard, Walter L., Jr.; Anderson, Melvin S.; Anderson, James Kent; and Card, 

Michael F.: Design, Analysis, and Tests of a Structural Prototype Viking Aero- 
shell. J. Spacecraft & Rockets, vol. 10, no. 1, Jan. 1973, pp. 56-65. 

13. Adelman, Howard M.; Lester, Harold C.; and Rogers, James L., Jr.: A Finite Element 

for Thermal Stress Analysis of Shells of Revolution. NASA TN D-7286, 1973. 


38 







(a) Geometry of cylindrical shell. 



(b) Nomenclature associated with finite-element 
representation of cylindrical shell. 

Figure 2.- Geometry and nomenclature of finite-element model. 


40 




(a) Cylindrical shell representation by 10 equally 
spaced elements. 


Region in 

Region n 
Region I 

Region I: ten elements - e/L = .025 
Region H: five elements - e/L » <,050 
Region in: five elements - e/L * .100 

(b) Conical shell representation by 20 unequally 
spaced elements. 

Figure 3.- Finite-element representations of shells. 
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Figure 4.- Comparison of methods at n = 0 and t = 1.2 see for a 
cylindrical shell having L/r = 2.4 and h/r = 0.05. 



rOIDIGNRL Q l STONGE > S/L 


Figure 5.- Effects of rotary inertia and transverse shear deformation on 
displacement at n = 0 and t = 1.2 sec for a cylindrical shell 
having L/r = 2.4 and h/r = 0.05. 
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— Linear solution 

-O Secondary effects 



Figure 6.- Effects of rotary inertia and transverse shear deformation on 
meridional moment resultant at n = 0 and t = 1.2 sec for a cylin- 
drical shell having L/r = 2.4 and h/r = 0.05. (MX = M x /(pv 0 ^r^).) 



Figure 7.- Effects of rotary inertia and transverse shear deformation on 
displacement at n = 0 and t = 1.2 sec for a cylindrical shell 
having L/r = 1.2 and h/r = 0.05. 



.651 



Figure 8.- Effects of rotary inertia and transverse shear deformation on 
displacement at n - 0 and t = 1.2 sec for a cylindrical shell having 
L/r = 0.3 and h/r = 0.05. 



Figure 9,- Comparison of higher order effects on displacement at n = 0 
and t = 2.4 sec for a cylindrical shell having L/r = 2.4 and 
h/r = 0.05. 
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MERIDIONAL DISTANCE. S/L 


Figure 10.- Effect of nonlinearities on displacement at n = 0 and 
t = 2.4 sec for a cylindrical shell having L/r = 2.4 and 
h/r = 0.025. 



MERIDIONAL DISTANCE. S/L 


Figure 11.- Effect of nonlinearities on displacement at n = 0 and 
t - 2.4 sec for a cylindrical shell having L/r = 2.4 and 
h/r = 0.0125. 
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Figure 15.- Mode shapes, modal stress resultants, and modal moment resultants for 

clamped-clamped cylinder. 



.7:383 
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(b) Second mode. 
Figure 15.- Continued. 
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(c) Third mode. 
Figure 15.- Continued. 
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